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Abstract 

Motivated  by  experimental  observations,  the  finite  element  method  is 
employed  to  model  the  competition  between  the  transverse  cracking  and 
delamination  modes  of  failure  that  occur  in  cross-ply  AS4/3501-6  gr/ep 
coupons  subjected  to  fatigue.  The  results  explain  the  extensive  delaminations 
and  reduced  crack  densities  that  arise  under  immersed  fatigue,  as  compared 
with  fatigue  in  air. 

1  Introduction 

The  effects  of  fluids  on  polymers  and  polymeric  composites  have  been 
studied  extensively  over  the  past  several  decades.  Recently,  the  subject 
received  a  new  impetus  due  to  the  revived  interest  in  utilizing  composites  in 
marine,  automotive,  and  civil  structures.  A  recent  review  of  the  literature 
(Weitsman,  1995)  indicates  that  relatively  little  information  is  available 
regarding  the  fatigue  response  of  composites  in  "wet"  conditions,  and  that  the 
sparse  amount  of  reported  data  concerns  pre-saturated  samples  fatigued  in  air. 

Recent  experiments  (Kosuri  and  Weitsman,  1995;  Gao  and  Weitsman, 
forthcoming)  on  the  fatigue  response  of  [0°/903]s  AS4/ 3501-6  gr/ep  coupons 
suggest  that  the  fatigue  life  of  pre-saturated  samples  subjected  to  cyclic  loading 
under  immersed  conditions  is  reduced  by  up  to  one  decade  in  comparison  with 
that  of  dry  coupons  fatigued  in  air.  Pre-saturated  coupons  fatigued  in  air  app>ear 
to  possess  the  longest  fatigue  life.  Perhaps  more  remarkably,  detailed 
observations  of  failed  coupons  have  shown  that,  when  compared  with  the  dry 
case,  the  saturated-immersed  circumstance  resulted  in  substantially  larger 
delaminations  at  the  0°/90°  interfaces  and  with  notably  fewer  transverse  cracks 
(i.e.  lesser  crack  density)  in  the  90°  ply  group.  Typical  observations  are  shown 
in  Fig.  1.  (Gao  and  Weitsman,  forthcoming). 

The  above  mentioned  experiments  utilized  simulated  sea  water.  They 
were  performed  at  a  frequency  of  5  Hz,  with  R  ratios  (R  =  amin/<ymax)  of  0.01,  0.1 
and  0.3  and  at  Omax  levels  of  400, 425,  450  and  475  MPa  (The  average  value  of  the 


coupons'  ultimate  failure  stress  was  540  MPa).  The  computations  below 
correspond  to  Omax  =  450  MPa  and  R=  0.1,  i.e.  Omm  =  45  MPa. 

Additional  experiments  were  performed  to  assess  the  rate  of  capillary 
climb  in  the  transverse  cracks  within  the  90°  ply  group  (work  by  Kosuri, 
reported  by  Smith  and  Weitsman,  1996).  Observations  suggest  that  such 

motion  proceeds  at  approximately  5  x  10'^  m/sec,  thus  requiring  about  250  sec.  to 

extend  over  the  entire  lengths  of  the  transverse  cracks  in  one  inch  (25  x  10'^  m) 
wide  coupons  immersed  in  sea  water.  This  time  span  exceeds  the  time-interval 
for  one  periodic  oscillation  in  the  fatigue  load  by  three  orders  of  magnitude. 

Thus  t  =  tc^illaiy/tfatigue  '*  10^- 

2  Basic  Considerations 

It  is  assumed  that  in  the  immersed  fatigue  case  water  enters  all 
transverse  cracks  by  capillary  motion  and  random  walk  process.  Subsequently, 
the  water  contained  within  the  cracks  cannot  escape  into  the  ambient  during 
the  down-loading  stages  of  the  fatigue  cycle  in  view  of  the  aforementioned 

value  of  t*=10^.  Thus,  the  water  trapped  within  the  transverse  cracks  is 
squeezed  into  all  weak  regions  within  the  composite  sample  diiring  the  down¬ 
loading  stages,  these  regions  being  the  0°/90°  interfaces.  Consequently,  the 
maximal  levels  of  delamination  fracture  energy  in  immersed  samples  occur 
when  the  pressure  p  on  the  crack  surfaces  is  at  its  peak  value,  namely  when 
^  =  Omin-  On  the  other  hand,  at  a  =  Omax  the  cracks  are  fully  open  and  p  =  o. 

The  above  reasoning  does  not  apply  to  fatigue  in  air,  in  which  case 
delaminations  are  attributable  to  stresses  near  or  at  a  =  Omax-  The  density  of 
transverse  cracks  is  associated  with  Omax  for  both  dry  and  immersed  cases, 
when  the  available  fracture  energy  for  that  failure  mode  is  largest. 

Turning  to  fracture  processes  it  is  hypothesized  that  all  transverse  cracks, 
which  span  the  inner  90°  ply-group,  are  equally  spaced  and  that  they  terminate 
just  short  of  the  0°/90°  interface.  It  is  also  assumed  that  all  delaminations 
emanate  from  the  tips  of  all  transverse  cracks  and  grow  in  unison  along  the 
0°/90°  interfaces,  as  shown  in  fig.  1. 

The  idealized  assumptions  of  equal  spacing  and  simultaneous  growth  of 
delaminations  are  made  to  retain  computational  tractability,  since  they  allow 
us  to  restrict  attention  to  a  p>eriodic,  "unit  cell"  geometry.  The  assumption  that 
transverse  cracks  do  not  impinge  on  the  0°/90°  interface  is  borne  out  by 
experimental  evidence.  This  assumption,  which  greatly  simplifies  the 
considerations  of  crack-tip  singularity,  can  be  further  justified  by  the  realization 
that  no  sharp  boundary  exists  at  the  0°/90°  interfaces  of  cross-ply  composites. 

In  addition,  it  is  considered  that  transverse  cracks  are  constituted  at  first 
across  the  thickness  of  the  90°  ply-group  and  subsequently  grow  in  a  straight 
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front  along  the  width  of  the  specimen  as  tunneling  cracks.  This  assumption  is 
also  borne  out  by  experimental  evidence. 

Finally,  the  analysis  is  restricted  to  the  evaluation  of  available  levels  of 
fracture  energies  in  several  static  states  and  does  not  consider  fatigue  effects  on 
crack  growth. 


3  Analysis 

3.1  Geometry  and  Material  Properties 


In  view  of  the  aforementioned  considerations,  attention  is  confined  to 
the  representative  cell  indicated  in  Fig.  2.  In  that  figure  h  =  0.5  mm  (half  the 
coupon’s  thickness)  and  L  =  12  mm.  The  length  L,  which  is  the  initial  value  of 
half  the  spacing  between  transverse  cracks,  was  chosen  after  some  preliminary 
computational  work  which  indicated  that  crack  interactions  are  insignificant  at 
spacings  that  exceed  L  =  12  mm.  The  coupon's  width  w=25mm  plays  no  role  in 
the  sequel,  since  we  consider  the  state  of  plane  strain.  Both  y  and  z  axes  in  Fig. 
2  serve  as  symmetry  lines  for  analysis  that  follows. 


The  Following  material  properties  and  hygro-thermal  coefficients  are 
assumed:  Ei  =  139GPa,  E2  =  E3  =1  IGpa,  V12  =  V13  =  0.27,  V23  =  0.34,  G12  =  4.8GPa, 
E3 


G23  = 


,  tti  =-  0.9|ie/“C,  a2  =  31m£/®C,  Pi  =  100pe/%Aw,  P2  =  3200|i£/%Aw. 


2(1+  V23) 

In  the  above  1  indicates  fiber  direction  and  fiber  direction  and  2  and  3  are 
mutually  perpendicular  to  it.  Ej,  Vy,  i,  j  =  1,2,3  are  the  stiffnesses  and  Poissons' 
ratios  respectively,  and  G12  and  G23  are  the  shear  moduli  in  the  1-2  and  2-3 

planes.  Ui’s  and  Pi’s  are  the  thermal  and  hygral  coefficients  of  expansion, 
respectively. 


Only  meager  information  is  available  regarding  the  effects  of  water  on 
the  levels  of  critical  fracture  energy  in  polymeric  composites  (Chou,  L,  et.  al., 
1993)  Consequently,  we  assume  that  the  same  values  apply  to  both  dry  and  wet 
cases. 


3.2  Finite  Element  Scheme 

All  computations  were  initiated  at  the  reference  configuration  shown  in 
Fig.  1.  Eight  node  quadrilateral,  plane  strain  elements  were  used  over  most  of 
the  domain  except  near  the  crack  tips,  where  six  node  triangular  plane  strain 
elements  were  incorporated.  PATRAN  [version  5,  MSC  1996]  was  used  to 
generate  the  mesh  and  develop  the  model  and  ABAQUS  [version  5.4,  Hibbit, 
Karlsson  and  Sorensen  Inc.  1994]  was  employed  to  solve  the  resulting  finite 
element  system. 


3 


As  noted  earlier  the  initial  spacing  of  transverse  cracks,  L  =  12  mm,  was 
considered  to  undergo  a  sequence  of  equal  subdivisions  under  fatigue  loading. 
At  each  stage  of  that  sequel  the  presence  of  delaminatidns,  as  well  as  their 
growth,  was  also  considered.  A  typical  circumstance  is  shown  in  Fig.  3, 
depicting  a  transverse  crack  spacing  of  2L  =  3  mm  and  delaminations  of  equal 
lengths  a.  In  accordance  with  the  basic  considerations  listed  earlier  the 

transverse  cracks  were  stopped  at  a  distance  5  =  7.8125  p  short  of  the  interface, 
which  corresponds  to  the  fiber  diameter.  Furthermore,  when  calculating 

energy  release  rates  along  the  0“/90°  interfaces  the  limiting  process  as  A  a  -4  0 

was  terminated  at  Aa  =0.18,  judging  this  to  be  the  smallest  reasonable  physical 
dimension  in  the  present  case.  Consequently,  the  smallest  elements  employed 

in  this  work  were  of  lengths  5,  and  0.15  as  noted  above. 

Finally,  it  is  worth  mentioning  that  in  all  circumstances  it  was 
ascertained  that  no  interpenetration  occured  across  opposing  crack  faces. 

3.3  Stress  States  and  Boundary  Conditions 

The  case  of  dry  fatigue  was  simulated  by  accounting  for  mechanical 
stresses  and  residual  thermal  stresses.  The  mechanical  stresses  were  evaluated 
under  peak  loads,  i.e.  Oy  =  Omax  =  450  MPa,  and  the  residual  thermal  stresses 

corresponded  to  a  temperature  excursion  AT  =  -125°C,  accounting  for  post¬ 
processing  cool  down.  In  that  case  the  surfaces  of  the  transverse  cracks  were 
considered  to  be  traction  free. 

Similar  considerations  applied  to  the  case  of  saturated  specimens 
fatigued  in  air,  except  that  the  residual  stresses  included  the  effects  of  a  tmiform 
moisture  profile  of  w  =  1.7%.  Note  that  all  surfaces  of  transverse  cracks  are  still 
taken  to  be  traction  free  and  the  largest  impetus  for  growth  of  both  damage 
mechanisms,  transverse  cracking  and  delaminations,  occurs  at  a  =  Omax- 

The  case  of  immersed  fatigue  differs  from  the  above  circumstances  in 
view  of  the  fact  that  the  water  entrapped  within  the  cracks  induces  a  hydrostatic 
pressure  on  all  crack  faces  at  all  levels  of  applied  stress  a  <  Omax*  This  pressure 
vanishes  at  a  =  Omax/  when  the  fully  open  cracks  can  freely  accommodate  the 
enclosed  water,  and  attains  its  largest  value  at  o  =  amn,.  Consequently,  while  the 
highest  impetus  for  further  transverse  cracking  still  occurs  at  a  =  aniax/ 
delamination  growth  is  driven  by  both  the  external  mechanical  stress  a  and  the 
now  present  internal  pressure  p.  Although  it  is  conceivable  that  maximal 
growth  could  occur  at  some  intermediate  level  of  a,  namely  aniin<o«ymax7 
corresponding  to  Pmax>P>o,  preliminary  computations  with  a  simpler  model 
(Smith  and  Weitsman,  1996)  suggest  that  it  is  reasonable  to  consider 

=  Omin  whereby  p  =pniax/  as  the  extreme  circumstance  for  the  instigation  of 
delamination  growth. 
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The  pressure  p  is  evaluated  by  the  iterative  procedure  outlined  below: 

As  a  first  step  the  volume  Vi  of  the  fully  open  cracks,  when 
o  =  o^max  andp  =  o,  is  evaluated  through  the  finite  element  analysis. 
Subsequently,  three  levels  of  hydrostatic  pressure  p  are  applied  to  the  crack 
surfaces  together  with  cj  =  <Tn,iii,  with  corresponding  crack  volumes  Vf 
computed  for  all  three  cases. 


The  results  are  correlated  by  a  quadratic  expression 

Vf  =  ap2  +  bp  +  c  (1) 


thus  a,  b,  and  c  are  know  constants. 


The  volumetric  change  between  Omax  and  Omm  is  related  to  p  through  the 
bulk  modulus  Kw  of  water  (Kw  =  2.5GPa),  thus 


P 


=  K 


‘W 


YLiYi 

Vi 


Expressions  (1)  and  (2)  yield  the  value  of  p  via 


(2) 


p2  +  p 


Kwb-Vj  c-Vj 
Kwa  a 


=  o 


(3) 


Obviously,  it  is  possible  to  determine  p  by  direct  iteration.  Comparative 
calculations  resulted  in  good  agreements  with  the  scheme  outlined  above, 
which  proved  to  be  much  more  efficient. 


4  Fracture  Energy  Analysis 

4.1  Transverse  Cracks 

As  noted  earlier,  these  cracks  are  assumed  to  grow  in  a  "tunneling" 
fashion  in  the  width  direction  of  the  coupon.  Considering  the  steady  state 
growth  of  new  cracks  mid-way  between  pre-existing  cracks  that  are  spaced  at 

distance  2X  apart,  straightforward  considerations  (e.g.  Weitsman  and  Zhu,  1993; 
Erdman  and  Weitsman,  forthcoming)  yield  the  following  expression  for  the 
available  fracture  energy  release  rate 

Q  _2U(X)-U(2X)  ,  . 

^  2c  ^  ^ 
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where  U(A,)  and  U(2A,)  denote  the  strain  energies  in  representative  cell  that 

correspond  to  crack  spadngs  of  X  and  2X,  respectively  and  2c  is  the  length  of  the 
crack  front,  (c  =  0.75h  in  the  present  analysis).  Note  that  only  mode  I  fracture 
energy  obtains  for  transverse  cracking. 

It  should  be  noted  that  the  crack  opening  displacements,  and  hence  also 
the  volume  enclosed  between  opposing  crack  faces,  may  diminish  upon  the 
formation  of  new  transverse  cracks.  This  implies  that,  for  expression  (4)  to 
hold,  some  of  the  water  entrapped  within  previously  established  cracks  is 
assumed  to  be  squeezed  out  during  ongoing  fatigue.  This  apparently 
reasonable  assumption  is  necessary  to  maintain  computational  tractability. 


4.2  Delaminations 

Delaminations  are  assumed  to  grow  in  unison  parallel  to  the  0“/90° 
interfaces  as  shown  in  Fig.  2.  They  are  considered  to  be  completely  established 
along  the  specimen’s  width. 

Mode  I  and  Mode  n  energy  release  rates,  and  are  evaluated 

employing  the  virtual  crack  extension  method  (Ribicki  and  Kanninen,  1977) 
and  its  subsequent  refinements  (Raju,  1987;  Raju,  Crews  and  Aminpour,  1988). 
The  above  quantities  are  computed  for  every  configuration  of  transverse  cracks 
and  for  delaminations  of  increasing  lengths  a  in  the  range  5  <  a  <  A,-5,  where  8 
denoted  the  fiber  diameter  as  before.  Note  that  the  presence  of  the  delimiting 

dimensional  scale  8  enables  us  to  disregard  the  oscillatory  singularity  at  the  tip 
of  the  delamination  and  thereby  determine  both  Gj  and  Gn- 

As  mentioned  earlier,  delaminations  in  the  case  of  immersed  fatigue 
were  evaluated  at  a  =  Omin,  which  corresponds  to  a  hydrostatic  pressure  p  =  pmax 
over  the  crack  surfaces.  This  circumstance  required  the  computation  of  Vi,  and 
the  subsequent  evaluation  of  p  =  Pmax  according  to  Eqs.  (l)-(3),  for  each 
configuration  of  transverse  cracks  and  at  all  delamination  lengths. 

4.3  Fracture  Criteria 

In  the  case  of  transverse  cracking,  new  cracks  are  assumed  to  occur  when 

Gi  =  Gic  (5) 

with  Gi  given  in  Eqn.  (4). 

Delamination  growth  is  considered  to  occur  when 
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Gdelam._  r'delam.  .  r^delam.  _  r* 

total  +^n  =^c  (o) 

where 

Gc  =  Gc(Gf G^^““  )  (7) 

as  suggested  recently.  (O'Brien,  1997). 

Obviously,  the  actual  values  of  and  G^*“"'  depend  on  the  spacing  of 

transverse  cracks,  as  well  as  on  the  delamination  length  a. 

The  competition  between  further  transverse  cracking  and  increasing 
delaminations  at  a  stabilized  configuration  of  transverse  cracks  is  determined 
by  comparing  criteria  (5)  and  (6).  For  prescribed  values  of  Oma,  and  Omin  the 
competition  settles  at  some  fixed  configuration  of  transverse  cracks,  beyond 

which  Gf®}™/Gc  >  Gi/Gic  and  delamination  becomes  the  preferred  mode  of 
failure. 

5  Results 

Results  for  the  dry  case,  and  for  the  fully  saturated  coupons  fatigued  in 
air  and  under  immersed  conditions,  are  exhibited  by  the  pairs  of  surfaces 
shown  in  Figs.  4,  5,  and  6,  respectively.  These  surfaces  correspond  to 
g  =  GavaUabie/Gcriticai  for  transverse  cracking  and  for  delaminations,  plotted  vs. 
the  density  of  transverse  cracks  and  of  delamination  length.  With  the 
exception  of  g  for  delamination  in  the  case  of  immersed  fatigue  (the  upper 
surface  in  Fig.  6),  which  corresponds  to  a  =  Omin  =  45  MPa,  all  other  surfaces 
comport  with  a  =  Omax  =  450  MPa.  Note  that,  as  a  consequence  of  our 
assumption  that  transverse  cracking  precedes  delamination,  the  progression  of 
transverse  cracking  is  independent  of  delamination  length  and  is  described  by 
the  developable  surfaces  in  Figs.  4, 5, 6. 

For  added  clarity  the  intersections  of  the  three  pairs  of  surfaces  in  the 
foregoing  figures  are  re-drawn  in  Fig.  7,  exhibiting  the  switching  from 
transverse  cracking  to  delamination.  Note  that  the  various  curves  in  Fig.  7 
emanate  from  disparate  levels  of  transverse  crack  densities  at  points  A,  B  and 
C,  namely  0.52  mm 0.89  mm-i  and  0.9  mm respectively.  In  all  circumstances 
the  initiation  of  delamination  is  followed  by  delamination  growth,  since 
G/Gc  >1  at  A,  B,  and  C.  However,  while  delamination  arrest  is  predicted  for 
fatigue  in  air  (curves  DIA  and  SI),  when  G/Gc  =  1,  unrestricted  growth  is 
predicted  for  the  saturated  immersed  case  (curve  SI).  The  calculated  density  of 
transverse  cracks  for  the  DIA  case  (0.89  mm-^)  compares  well  with  experimental 
observations  (1.0  mm-^).  However,  experimental  observations  indicate  that  a 
higher  crack  density  (0.3  mm-i)  occurs  in  the  SI  case  than  predicted  in  the 
present  article  (i.e.  0.52  mm-i).  The  disparity  suggests  that  Gc  in  the  immersed 
case  is  lower  than  that  in  the  dry  case,  in  contrast  to  the  presumed  equality 
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employed  in  the  calculations.  However,  no  Gc  data  are  available  for  the 
immersed  case. 

6  Concluding  Remarks 

A  detailed  finite  element  analysis  was  performed  to  evaluate  the 
competition  between  the  progression  of  transverse  cracking  and  interfacial 
delaminations  in  cross-ply  gr/ep  coupons  subjected  to  fatigue.  The  analysis 
considered  the  circumstances  of  ktigue  in  air,  for  both  dry  and  water-saturated 
samples,  and  of  the  immersed  fatigue  of  saturated  coupons.  An  accounting  for 
fatigue  was  attained  by  evaluating  available  fracture  energies  at  the  maximal 
and  minimal  load  levels  during  cycling. 

To  maintain  tractability  the  analysis  employed  idealized  geometric 
configurations  that  assume  symmetric  formations  of  both  transverse  cracks  and 
delaminations.  While  such  idealizations  tend  to  overestimate  all  levels  of 
available  fracture  energy,  when  compared  with  realistic  circumstance  which 
involve  fewer  cracks,  this  trend  may  be  compensated  by  the  fact  that  the 
bending  which  arises  in  the  absence  of  symmetry  may  greatly  enhance  the 
foregoing  energy  levels. 

It  is  interesting  to  note  that  the  model  and  results  presented  herein, 
which  involved  static  analyses,  agree,  at  least  qualitatively,  with  experimental 
observations. 
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Fie.  1.  Observations  of  Transverse  cracks  in  the  90°  Ply-Group  Delaminations  at  the  0°/90°  Interfaces  of  Failed  [0°/903]s  gr/ep 
Coupons  Subjected  to  Fatigue.  The  Peeled-Off  Portion  of  the  0°  Ply  is  shown  on  the  Top  of  Each  Photo,  Indicating  Extent 
of  Delamination,  (a)  Dry  and  (b)  Saturated  Immersed  for  amax  =  0.8au.  (c)  Dry  and  (d)  Saturated  Immersed  for 

Omax  =  O.SSOu.  In  all  cases  R  “  ^min/^max  —  0.1. 


Fig.4.  Available  Fracture  Energies,  scaled  by  the  required  levels,  for 
Transverse  Cracking  and  for  Delamination  at  the  maximal  stress  level 
(a=450MPa)  of  Dry  Fatigue,  vs.  Density  of  Transverse  Cracks  and 
Delamination  Lengths  in  a  [0®/903]s  gr/ep  coupon. 
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Fig  5.  Available  Fracture  Energies,  scaled  by  the  required  levels,  for 
Transverse  Cracking  and  for  Delamination  at  the  maximal  stress  level 
(a=450MPa)  of  saturated  coupons  fatigued  in  air,  vs.  Density  of 
Transverse  Cracks  and  Delamination  Lengths  in  a  [0°/90§]s  gr/ep 
coupon. 


Fig.  6.  Available  Fracture  Energies,  scaled  by  required  levels,  for  Transverse 
Cracking  at  the  maximal  stress  level  (o=450  MPa)  and  for  Delamination 
at  the  Downloading  stage  (<j=45  MPa)  of  Immersed  Fatigue  vs.  Density 

of  Transverse  Cracks  and  Delamination  Lengths  in  [07903]s  gr/ep 
coupons. 


CRITICAL  DENSITY  OF  TRANSVERSE  CRACKS 
(A)  0.52  mm'^  (B)  0.89  mm'^  (C)  0.9  mm"^ 


Fig.  7.  The  Intersection  Curves  of  the  surfaces  shown  in  Figs.  5,  6,  and  7, 
indicated  by  DIA  ("Dry  in  Air"),  SIA  ("Saturated  in  Air"),  and  SI 
("Saturated  immersed"),  respectively.  These  Curves  exhibit  the  trends 
towards  increasing  delaminations  upon  switching  from  failure  by 
transverse  cracking.  Note  that  the  above  curves  emanate  from  distinct 
values  A,B,C  of  transverse  crack  density.  Also,  Delamination  Arrest  is 

predicted  for  the  DIA  and  SIA  cases,  when  G/Gc^l,  while  no  arrest  is 
predicated  for  the  SI  case. 


